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ABSTRACT 

The Long Duration Exposure Facility (LDEF) has provided an 
excellent opportunity to understand the nature of directed atomic 
oxygen interactions with protected polymers and composites. 
Although there were relatively few samples of materials with 
protective coatings on their external surfaces on LDEF which were 
exposed to a high atomic oxygen fluence, analysis of such samples 
has . enabled an examination of the shape of atomic oxygen undercut 
cavities at defect sites in the protective coatings. 

Samples of front-surface aluminized (Kapton) polyimide were 
inspected by scanning electron microscopy to identify and measure 
crack defects in the aluminum protective coatings. After 
chemical removal of the aluminum coating, measurements were also 
made of the width of the oxidized undercut cavities below the 
crack defects. The LDEF flight undercut cavity geometries were 
then compared with Monte Carlo computational model undercut 
cavity predictions. 

The comparison of the LDEF results and computational modeling 
indicates agreement in specific undercut cavity geometries for 
atomic oxygen reaction probabilities dependent upon the 0.68 to 
3.0 power of the energy. However, no single energy dependency 
was adequate to replicate flight results over a variety of 
aluminum crack widths. 


737 


INTRODUCTION 


Atomic oxygen in low Earth orbit readily oxidizes most 

polymericmaterials ( 1-5 ) . 

It is now widely known that the use of unprotected 

polymericmaterials in low Earth orbit may ca " s ®. ^ CCe ^^eJ new 
recessionduring high atomic oxygen fluence missions. Either ne 
materialsshould be used which are inherently immune to atomic 
oxvgenattack or conventional materials may continue to be used 
providedthey are covered with atomic oxygen resistant coatings. 

Protective coatings for polymers typically used in low 
Earthorbit, such as SiO x on polyimide Kapton, have been developed 
andextensively tested in ground laboratory atomic oxygen systems. 
However, few high fluence in-space exposure studies of 
suchprotected polymers have been conducted ( 7 , a ) . 

On the LDEF spacecraft, relatively few polymers that had metal 
ormetal oxide protective coatings on their exposed surf a 
wereexposed to high fluence atomic oxygen. A small sample 
of aluminized (Kapton) polyimide was exposed to high fluence 
atomicoxygen on LDEF and then used for characterization of th 
undercutcavities at defect sites in the protective coating 
These datawere then used for comparison with an ab initio Monte 
SriocompStStional model which was developed to predict the shape 
ofthe atomic oxygen undercut cavities that occur at defects 
innrotective coatings on polyimide Kapton. The 

detailedmechanistic assumptions of the L DEF 

then bevaried to examine the degree to which the observed LDEF 
resultsmatched the model predictions for various 
mechanisticassumptions . 

Understanding the detailed mechanistic interactions which 
occurbetween atomic oxygen and protected materials in low Earth 
orbitwill allow more reliable predictions °f the *° ng 
durabilityof protected materials in low Earth orbit, based on 

groundlaboratory testing . 

METHODS, MATERIALS, AND PROCEDURE 

T.DEF Aluminized Kapton Samples 

A front-surface aluminized (Kapton H) polyimide sample exposed 
onLDEF as part of the Solar Array - Materials Passive 
LDEFExperiment located on row 8A was used for the experimental 
datapresented in this investigation. The sample consisted 
of approximately 1,000 k of vapor deposited aluminum on J CK0254 
mmthick Kapton substrate. The sample was oriented such that 
itreceived a total atomic oxygen fluence of 7.15 x lu 
atoms /cm 2 arriving at 38.1° from the perpendicular to the surface. 
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The sample was examined for crack or scratch defects in the 
aluminization to allow a two-dimensional characterization of the 
crack width and the associated undercutting. A two-dimensional 
geometry was desired because the Monte Carlo computational model 
was also two-dimensional. 

Scanning electron microscope photographs were taken of defects 
with a variety of widths in the aluminized Kapton. The sample 
was then chemically treated with a dilute hydrochloric acid 
solution to remove the aluminum. Additional scanning electron 
micrographs were taken at the identical defect sites to allow 
measurement of the width of the undercut cavities associated with 
the defect sites. 

The crack-width data, undercut-width data, and undercut-cavity 
shape were then used as the basis for modeling this behavior with 
Monte Carlo predicted undercutting geometries based on specific 
atomic oxygen interaction mechanistic assumptions such as 
reaction probabilities and fractional energy loss upon impact. 

Monte Carlo Computational Model 

An ab initio Monte Carlo computational model was developed to 
predict the undercut cavity shape as a result of atomic oxygen 
interaction with polyimide Kapton at the site of crack defects in 
the protective aluminum coating. The model allows the atomic 
oxygen to enter the crack from a fixed direction (as was the case 
for LDEF ) , a sweeping direction (as would occur in the case of a 
solar array surface where the surface rolls with respect to the 
atomic oxygen arrival direction), or an isotropic arrival 
direction (as occurs in ground laboratory oxygen plasma exposure 
tests to simulate low Earth orbital atomic oxygen conditions). 

For the fixed direction arrival case, the model also allows a 
thermally distributed arrival direction to be assumed because the 
atoms are hot (=1227K for the LDEF mission) (8), which causes 
transverse velocity variations to occur. 

The mathematical model of Kapton underneath the crack defect 
consisted of an orthogonal array of cells in which the details of 
the atomic oxygen interaction with each cell were based on 
mechanisms which were mathematically defined. Figure 1 
illustrates the geometry of the Monte Carlo computational model. 
Figure 2 is a flow chart of the Monte Carlo computational model 
showing the decision-making process and the iterations used to 
predict the undercut cavity shape. 

The location of atoms being brought into the defect in the 
aluminum coating is randomly selected. When the atomic oxygen 
model atom impinges upon the Kapton surface, a random selection 
is made as to whether or not the oxygen reacts with the Kapton. 
This selection is based on a prescribed initial impact reaction 
probability. In general this reaction probability was set at 
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0.138 (for normal incidence) based on previous experimental and 
computational investigations ( 9 ) . 

The probability of reaction upon energetic impact was also 
assumed to be dependent upon the arrival angle, 0, with respect 
to the perpendicular to the surface with a (cosine @) 1/2 
dependence (9,10). The angle-dependent reaction probability was 
based on the impinging direction with respect to the local 
surface normal which was based on the average of adjoining Monte 
Carlo computational model cells. 

Oxygen atoms which impinged upon aluminum were treated as if a 
glass surface was struck where no reaction was possible. However, 
depending upon the program operator's choice, recombination or 
ejection and scattering distributions could be prescribed as 
similar to those used for atomic oxygen ejection from Kapton. 

Energy loss upon impact was defined in terms of a fraction of 
energy loss between the arriving energy and that energy which 
would be associated with thermally accommodated atomic oxygen. 
Thus, the atomic oxygen was assumed to be at 4.5 eV upon initial 
impact and ultimately thermally accommodated down to 
approximately 0.04 eV after numerous collisions. The atomic 
oxygen reaction probability was explored for a dependence upon 
energy proportional to the arriving energy raised to a numerical 
exponent ranging from 0.68 to 4. 

The atomic oxygen which did not react was ejected at an angle 
which could be selected based on choices made by the Monte Carlo 
program operator. These choices included elastic scattering, 
cosine distributed scattering, or a combination of elastic 
scattering and thermally accommodated cosine distributed 
emission. One could also select a random probability of 
recombination of atomic oxygen if it did not react with the 
Kapton . 

The Monte Carlo computational model ray traced each atom until it 
either reacted with Kapton or exited the undercut cavity. 

Exiting could be accomplished by leaving through the defect in 
the aluminum, or, in the case of high fluence exposure, exiting 
from an aperture produced in the bottom of the Kapton which was 
not protected. 

Undercutting profiles were predicted by the Monte Carlo model and 
measurements were then made from the undercut profiles produced 
by the model and compared with actual LDEF data. 

The use of the Monte Carlo computational simulation requires a 
calibration in which actual in-space erosion for a given atomic 
oxygen fluence is correlated with computational modeling erosion. 
This is accomplished by using a wide defect in the protective 
coating to simulate the erosion which would occur in space, using 
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the Monte Carlo computational model. 

The mean depth of erosion for simulated unprotected Kapton was 
then used to predict the number of atoms, N, required to simulate 
the actual LDEF in-space damage which would occur at the site of 
a narrow crack for which experimental data exists from LDEF. 

Thus the number of atoms required to enter the simulated defect 
to cause the same damage as would occur during the LDEF mission 
is given by 


where : 


N = MEFW 
HDC 


M = number of atoms 

entering a wide defect for computation simulation 
calibration 

E = in-space erosion yield of polyimide Kapton, cmVatom 

F = atomic oxygen fluence in space for the LDEF coupon, 
atoms /cm 2 

W = narrow crack width in Monte Carlo cell units representing 
the LDEF protective coating cracks 

H = wide crack width in units of Monte Carlo cells for 

simulation of unprotected Kapton for Monte Carlo calibration 

D = depth of erosion in units of Monte Carlo cells for the wide 
defect calibration Monte Carlo simulation 


C = Monte Carlo simulation cell edge length, cm/cell. 

With the above equation one could then use the Monte Carlo model 
to predict the atomic oxygen undercutting which would occur in 
space. The above equation also allows calculation of the scale 
for the Monte Carlo computational geometry. In the above 
equation a Kapton in-space erosion yield, E, of 3 x 10' 24 cmVatom 
was assumed. 


RESULTS AND DISCUSSION 

The aluminized Kapton sample retrieved from LDEF had cracks which 
were oriented in such a direction that if one looks down the 
length of the crack, the average angle of atomic oxygen arrival 
was 25.9° from the perpendicular to the coated Kapton surface. 
Figure 3 shows scanning electron micrographs of this sample at a 
crack defect site in the aluminized coating, before and after 
removal of the aluminized coating. 
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As can be seen from Figure 3, the width of undercutting 
significantly exceeds the width of the crack defect in the 
aluminized coating. 

Figure 4 shows a clear view of the undercut cavity of a similar 
defect site after removal of the aluminized coating. As can be 
seen from Figure 4, the undercut cavity is asymmetrical. The 
left side of the undercut cavity has a chamfer near the top while 
the right side of the undercut cavity, which does not contain a 
chamfer, is undercut at such an angle that the cavity wall is not 
visible in this image. This is thought to be largely due to the 
fact that incoming atomic oxygen scatters off the edges of the 
aluminized coating causing the scattered atomic oxygen to 
preferentially attack the near surface Kapton opposite the side 
of the defect. 

In Figure 4, the atomic oxygen was arriving from the upper left 
direction. Thus the undercut cavity has rather straight walls on 
its right side and is more undercut on the left side because of 
atomic oxygen scattering effects. It was apparent from the 
micrographs of wide cracks (of the order of 3^m on this sample) 
that the atomic oxygen undercut cavity extended all the way 
through the 0.0254 mm thick Kapton sample. Narrow cracks in the 
aluminized layer did not allow sufficient atomic oxygen fluence 
entry into the undercut cavities to erode all the way through the 
Kapton, as shown in Figure 5. 

Figure 5 also shows the chamfering effects of atomic oxygen 
reflected from the aluminized coating. The bottom of the 
undercut cavity is also observable at this rather narrow crack 
defect site. 

Figure 6 shows the results of measurements from the scanning 
electron micrographs of a variety of crack widths in the 
aluminized coating and their associated undercut widths , for . 
cracks all oriented in approximately the same direction. This 
plot, along with the asymmetrical shape of the undercut cavities 
observed in the scanning electron micrographs, was used as a 
basis of comparison with Monte Carlo model predictions to assess 
the validity of mechanistic assumptions for atomic oxygen 
interaction with polymers and protective coatings. 

The Monte Carlo modeling of the LDEF atomic oxygen arrival used 
the assumptions of reference (8). The angular distribution of 
atomic oxygen arrival associated with 1227K hyperthermal atoms is 
shown in Figure 7 as a polar plot of the atomic oxygen arrival 
flux as a function of the arriving direction. 

As can be seen from Figure 7, the arriving atomic oxygen has an 
angular distribution similar to a searchlight which clearly would 
indicate that undercut cavities should widen with depth, if a 
sufficient atomic oxygen fluence has passed through the defect. 
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Although computational simulations of the LDEF environmental 
interactions with aluminized Kapton were conducted assuming a 
variety of atomic oxygen scattering processes (including 
specular, cosine, and a mix of cosine and specular), only a 
cosine ejection distribution was found to produce undercut 
profiles which replicated observed LDEF results. 

Figure 8 shows two Monte Carlo predicted undercut cavities. 

These two predictions were selected for providing reasonable 
undercut cavity matches to both the undercut shape and the 
undercut cavity width versus defect width plot of Figure 6. Such 
models were then developed for reaction probability energy 
dependent exponents ranging from 0.68 to 3. 

As can be seen from Figure 8, atomic oxygen reflected from the 
aluminized protective coating does produce a slight pocket in the 
Kapton near the aluminum coated surface as shown in the observed 
LDEF scanning electron micrographs. 

For each reaction probability energy-dependent exponent, a 
variety of fractional energy losses upon impact were evaluated in 
an attempt to reproduce the LDEF-observed undercut cavity width 
versus protective coating defect width data. The results of 
these efforts are shown in Figure 9. The Monte Carlo predictions 
are shown as error bars on the figure. The error bars represent 
the range of undercut widths that could be predicted based on the 
range of interpretations possible for the width of the 
undercutting profile. No individual reaction probability energy 
dependent exponent and its associated fractional energy loss upon 
impact was found to produce an ideal match to the observed 
experimental data. 

A rather weak energy dependence as shown in Figure 9a produces 
far too much undercutting for small width defects. Increasing 
the energy dependence up to an exponent of 3, Figure 9d did 
produce crack width to undercut width ratios which replicated 
flight results; however, predicted undercutting for large width 
cracks was not sufficiently wide to match observed data. Thus, 
although undercut widths could be predicted which matched LDEF 
results for specific width cracks, no single energy dependent 
exponent for reaction probability was found to provide acceptable 
profile matches over all width cracks. Based on these results 
further adjustments in the quantification of the mechanistic 
assumption paraments are necessary. 

One of the difficulties associated with performing the Monte 
Carlo computational calculations is that typically ten hours were 
needed to produce each predicted Monte Carlo undercutting 
profile. This is largely a result of Monte Carlo model atoms 
bouncing around within the undercut defect cavities with a very 
low reaction probability after they have thermally accommodated. 
Each atom is required to either react or exit the defect before 
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another atom is brought into the model cavity. No recombination 
was assumed to occur in these studies. As the reaction 
probability decreases with each collision, many impacts are 
required before the next atom is allowed to enter. Thus, the 
process of computation slows down greatly as the undercut cavity 
grows. Although this is an inconvenience, the accuracy of the 
prediction should not be compromised. 

The Monte Carlo model does predict many of the features 
experimentally observed but will require additional investigation 
to more clearly quantify the energy dependence and fractional 
energy loss. 


SUMMARY 


A retrieved LDEF sample of front surface aluminized Kapton was 
analyzed by scanning electron microscopy to characterize the 
atomic oxygen undercut geometries and the dependence of undercut 
cavity width on the aluminized protective coating crack width. 

The undercut cavities observed indicate that atomic oxygen 
scatters off the edges of the protective coating to produce 
accentuated atomic oxygen erosion of the Kapton opposite the 
impinged edge. A Monte Carlo computational model was developed 
which is capable of simulating a variety of atomic oxygen 
environments including that of the LDEF spacecraft. This model 
was used to predict undercut cavity profiles in an attempt to 
replicate the observed LDEF data. 

Although the results of the computational modeling were found to 
replicate the shape of some of the undercut cavities at defect 
sites, such replication was not observed over a wide range of 
crack defect widths. Further study is needed to develop 
mechanistic information which accurately predicts observed 
experimental results . 
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Figure 1. 


Monte Carlo computational model geometry . 




GENERAL OUTLINE 



Figure 2 . 


Monte Carlo computational model flow chart. 











a. A1 coating present 



b. A1 coating removed 


Figure 3. Scanning electron photomicrograph of retrieved LDEF 
aluminized Kapton sample 
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Figure 4. Undercut cavity associated with a wide defect in the 
aluminized Kapton after removal of the aluminized coating. 



Figure 5. Undercut cavity beneath a narrow crack in the 
aluminized coating after removal of the aluminized coating. 
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15.0 



Figure 6. Observed retrieved LDEF sample undercut cavity width 
dependence upon crack width in the aluminized coating 



4 

SPACECRAFT RAM DIRECTION 


Fiaure 7. Polar plot of LDEF atomic oxygen arrival resulting 
from a combination of the spacecraft's orbital velocity. Earth's 
upper atmosphere co-rotational velocity, the 28.5 inclination of 
the orbit, and the fact that the spacecraft is ramming into a 
1227K gas of atomic oxygen. 
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a * ^’o^68 lin w ^ e cra ck defect, reaction probability dependent upon 
(energy) , fractional energy loss upon impact equal to 0.4. 



b. 2 urn wide defect in the aluminized coating, reaction 
probability dependent upon (energy) 3 , fractional energy loss upon 
impact equal to 0 . 2 . 


^*^-?fure 8. Monte Carlo computational model predictions for undercut 
cavities associated with crack defects in LDEF aluminized Kapton. 
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Undercut Width (microns) Undercut Width (microns) 


Crack Width (microns) 

a. Reaction probability proportional to (energy) 
:tional energy loss upon impact = 0.4 


Crack Width (microns) 


b. Reaction probability proportional to (energy) 1 , 
fractional energy loss upon impact = 0.2 


Figure 9. Comparison of Monte Carlo computational model 
predicted undercutting (shown as error bars) with experimentally 
observed LDEF data. 
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c. Reaction probability proportional to (energy) 2 , 
fractional energy loss upon impact =0.2 



d. Reaction probability proportional to (energy) 3 , 
fractional energy loss upon impact =0.2. 


Figure 9. Comparison of Monte Carlo computational model predicted 
undercutting (shown as error bars) with experimentally observed LDEF 
data. 
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